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I. INTRODUCTION 



Nuclear structure changes significantly when ap- 
proaching closed shells. The change of the nuclear shape 
from deformed to spherical is expected to influence the 
level density as well as the 7-strength function, which 
are important in order to calculate cross-sections ^ and 
(n,7) reaction rates |^]. These quantities are critical in- 
put parameters in astrophysical models describing the 
nucleosynthesis in stars 

The Oslo Cyclotron group has developed a method to 
simultaneously extract the level density and 7-strength 
function, over a wide excitation energy and 7-energy re- 
gion, from measured 7-ray spectra |^,^. Today in the 
rare earth region only mid-shell nuclei have been investi- 
gated ■ Therefore a natural extension of these stud- 
ies is to pin down the properties for nuclei close to the 
TV = 82 gap. This work presents experimental data on 
level densities and 7-strength functions of the weakly 
deformed ^'^^Sm and ^^^Sm nuclei. The present exper- 
imental technique represents the least model-dependent 
method to obtain the 7 strength function over a wide 
energy region below the neutron-separation energy. It 
is actually the density of levels accessible to the nuclear 
system, in the 7-decay process, which is extracted ex- 
perimentally and interpreted as the level density. This 
interpretation is approximately valid for a thermalized 
nucleus in the continuum. The main advantage of utiliz- 
ing 7-rays as a probe for level density is that the nuclear 
system is likely to be thermalized prior to the 7-emission. 

A long-standing problem in experimental nuclear 
physics has been to observe the transition from a super- 



fluid state, at the ground state, to a normal (Fermi gas) 
state at higher temperature. The nuclear level density 
can be utilized to deduce e.g. entropy, temperature and 
heat capacity and these observables can reveal indica- 
tions of a phase transition. 

The 7-strength function is a measure for the average 
electromagnetic properties of nuclei and has a fundamen- 
tal importance for the understanding of nuclear struc- 
ture and reactions involving 7 rays @. At - 3 MeV 
of 7 energy, a bump is observed in the 7-strength func- 
tion of rare earth nuclei from (^IIe,a) experiments [^j. 
The structure in the 7-ray strength function, usually 
called pygmy resonance, will be discussed in terms of the 
scissors mode built on exited states. The samarium re- 
sults are compared to results for well deformed ^^'^•^^'^Dy^ 
166,167jt;j. g^^^ i7i,i72yi-, isotopes, with data from (n,7)- 
experiments |lll and giant dipole resonance data. 

Section II describes the experimental methods. In 
Sect. Ill the experimental level density and 7-strength 
function of ^**Sm and ^"^^Sm are obtained. Section IV 
examines thermodynamic properties within the micro- 
canonical and the canonical ensemble, while electromag- 
netic properties of the two nuclei are discussed in Sect. V. 
Conclusions are made in Sect. VI. 



II. METHOD 



The experiment was carried out at the Oslo Cyclotron 
Laboratory (OCL). The neutron pick-up (^He,a) and the 
inelastic scattering (•^IIe,^He')-reactions where employed 
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on a -'^^^Sm self-supporting target, with the 45 MeV ^He- 
beam dehvered by the MC-35 cyclotron. In the inelas- 
tic scattering reaction collective excitations are probably 
more pronounced than in the pick-up reaction, however, 
in Ref. it is shown that the two reactions give similar 
results. 

Charged particles and 7-rays were recorded with the 
multi-detector system CACTUS , which contains 8 Si 
particle telescopes and 27 Nal 7-ray detectors. The par- 
ticle telescopes are placed at an angle of 45° relative to 
the beam axis, and consist of a Si front (AE) and a Si(Li) 
back (E) detector with thickness 140 and 3000 /im, re- 
spectively. The Nal 7-detector array, having a resolution 
of ~ 6 % at 7 energy = 1 MeV and a total efficiency 
of ~ 15 %, surrounds the target and particle detectors. 
In addition, two Ge detectors were used to monitor the 
spin distribution and selectivity of the reactions. The 
typical spin window in the two reactions employed is 2-6 

h. 

From the reaction kinematics the measured ejectile en- 
ergy can be transformed to excitation energy E. Thus, 
each coincident 7-ray can be assigned to a 7-cascade orig- 
inating from a specific excitation energy. The data are 
sorted into a matrix of {E^E^) energy pairs. Examples 
of the recorded 7 spectra, the so called raw spectra, from 
5 MeV excitation energy are shown in the left panels of 
Fig. for "^Sm (top) and ^^^Sm (bottom). Note that the 
statistics of ^''^Sm is about twice as good as for ^'"^Sm, 
since the (^IIe,'^IIe') reaction has a higher cross section 
than the ('^IIe,a) reaction. For each excitation energy E 
the Nal 7-ray spectra are corrected for the measured re- 
sponse function of the Nal detectors with the unfolding 
procedure of Ref. |Q. Unfolded 7 spectra for the two 
nuclei are shown in the middle panels of Fig. |l|. 

From the corrected 7-ray spectra the primary 7-ray 
matrix was extracted according to the subtraction tech- 
nique ofRefs. @,|l|l. Examples of primary 7 spectra can 
be seen in the right panels of Fig. |l|. The method of ex- 
tracting the primary 7 spectra is based on the assump- 
tion that the decay properties of a bin of excited states 
are independent of whether the states are directly popu- 
lated through the nuclear reaction or from de-excitation 
from higher excited states. This is believed to be ap- 
proximately fulfilled because the long life time of excited 
states gives the nucleus time to thermalize prior to the 7 
decay. 



III. EXTRACTION OF LEVEL DENSITY AND 
7-STRENGTH FUNCTIONS 

The energy distribution of primary 7-rays provides 
information on both the level density and the 7-ray 
strength function, enabling a simultaneous determina- 
tion of the two functions. The fundamental assumption 
behind the extraction procedure is the Brink-Axel hy- 
pothesis [p^|Jl7|| , where the probability of 7-decay in the 



statistical regime, represented by the primary 7 matrix 
P{E,E^), can be expressed simply as a product of the 
final-state level density p{E — Ej) and a 7-energy depen- 
dent factor F{E-y) 



P{E, E^) cx F{E^)p{E - E^). 



(1) 



From the 7-energy dependent factor F{E^) the 7- 
strength function can be found. The details of the 
method and the assumptions behind the factorization of 
this expression are described in Ref. |^ and only an out- 
line of the normalization procedures is given here. Equa- 
tion (|l|) has an infinite number of solutions. It can be 
shown that all equally good solutions of Eq. ([l]) can 
be obtained by the transformations 

p{E - E^) = Aexp[a(£; - E^)]p{E - E^), (2) 
F{E-y) = Beyii>{aE-y)F{E-y) (3) 

of any particular solution (p,F). Consequently, neither 
the slope nor the absolute value of the two functions 
can be obtained through the iteration procedure, but the 
three variables A, B and a of Eqs. (|^) and (H) have to 
be determined independently to give the best physical 
solution of the level density and 7-strength function. 



A. Level density 

The parameters A and a of Eq. (|^) can be determined 
by fitting the level density from the iteration procedure 
Q to the number of known discrete levels at low ex- 
citation energy and to the level density estimated from 
neutron-resonance spacing data at high excitation en- 
ergy |19|. This normalization procedure is shown for 
^^^Sm in Fig. |^. The deduced level density is fitted to the 
discrete levels (between the arrows in the upper panel of 
Fig. HI) as far up in energy as we can assume that all levels 
are known. At high excitation energies the deduced level 
density is fitted (between the arrows in the lower panel 
of Fig. |2|) to a Fermi-gas approximation of the level den- 
sity (line). The Fermi-gas approximation is forced (multi- 
plied by a constant / given in Table |) to pass through 
the level-density estimate at the neutron-separation en- 
ergy (filled square) obtained from the neutron-resonance 
spacing data. The level density at the neutron binding 
energy Sn is estimated using the Fermi-gas expression: 



PSn 



D (/ + 1) exp(-(J-M)2/2a2) + /exp(-/V2cr2) 

(4) 



where D is the neutron resonance spacing, / is the spin 
of the target nucleus in neutron resonance experiments 
(^i^'SSm) = a and I("^Sm) = 7/2h), the spin-cuto ff pa- 
rameter a is defined by a'^ ^ 0.0888 v/a(5„ - Ets)A'^/^, 
while A is the mass number, and Eb^ the back-shift pa- 
rameter p^ ]. Some of the parameters used are given in 
Table |. 
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The level densities for -'^'^^Sm and ^''^Sm as a function 
of excitation energy are shown in Fig. ^. The fact that 
the even nuclei has lower level density is well understood 
as due to pairing making the even nuclei more bound by 
bringing the ground state down in energy. The physical 
interpretation of these data is given in Sect. IV. 



B. Gamma-strength function 

Blatt and Weisskopf |Q suggested to adopt the ratio 
of the partial radiative width Ti{E^) and the level spac- 
ing of the initial states i with equal spin and parity Di 
in order to describe the 7 decay in continuum. Accord- 
ing to this, the definition of the 7-ray strength function 
IS given by /xL = Ti{E-,)/{E^^+^D,), where X denotes 
the electric or magnetic character, and L defines the mul- 
tipolarity of the 7 transition. After the normalization of 
the level density, the parameter a of Eq. (||) is already 
fixed, while the normalization constant B remains to be 
determined. The 7-energy dependent factor F{E^) is 
proportional to 



E 

XL 



(5) 



where /xL^Ey) is the 7-ray strength function for the 
multipolarity XL. We assume that the 7 decay in the 
continuum of nuclei with low spin is mainly governed by 
electric and magnetic dipole radiation and that the acces- 
sible levels have equal numbers of positive and negative 
parity states. Thus, the observed F{E^) can be expressed 
by a sum of the E\ and M 1 7-strength functions only: 



BF{E-,) = [fEiiE-,) 



fMi{Ej)]E^. 



(6) 



The average total radiative width of neutron resonances 
(r^(S'„, /, tt)) with excitation energy equal to the 
neutron-separation energy 5„, spin /, and parity n 



■EE 



p{S„,I,7r)——^j, 



dE^E^^+^ 



fxL{E^)p{Sn - E^,If,TTf), 



(7) 



can be written in terms of F{E^) by means of Eq. (|^) 
for L — I. The experimental value of the average to- 
tal radiative width (F^) is the weighted sum of Eq. (0) 
with contributions from all different / accessible in the 
(n, 7) reaction. Thus, with the experimental level density 
p already normalized, the normalization constant B can 
be deduced. The average total radiative width of neu- 
tron resonances used for the normalization of the ^'^^Sm 
7-strength function is (F^) = 69(2) meV, as given in 
Ref. ^ and for "^Sm (F^) = 45 meV from Ref. Q is 
used. A detailed description of the normalization of the 
7-strength functions is given in Ref. S . The normalized 
7-strength functions of ^"^^Sm and ^ Sm are shown in 



Fig. IJ. The observed strength functions differ from each 
other with less than a factor of two. 

In order to check the quality of the level density and 
7-strength functions, the total 7-spectra for ^'^^Sm and 
^■^^Sm (originating from excitation energy equal to the 
neutron binding energy) has been calculated and com- 
pared with experimental ones in Fig.jsj. The data points 
with error bars are taken from the ^'^^m(n,7)^'*®Sm and 
^^^Sm(n,7)^''^Sm experiments The solid line is cal- 
culated from the 7-strength functions and level densi- 
ties extracted from the present ^'**Sm(''lIc,a7)^^^Sm and 
i**^Sm(3He,3He'7)"*Sm data. The calculation is per- 
formed by averaging over 100 keV intervals. The fact 
that the agreement it so good gives added confidence that 
our method is working well. 



IV. THERMODYNAMICAL NUCLEAR 
PROPERTIES 

The present experimental technique has been previ- 
ously applied to several well deformed rare earth nuclei. 
In Fig. ^the level density functions for ^''^'^''^Sm are com- 
pared with the level density functions for the I6i,i62py 
||, 166467EJ. Q 171472 Yb g nuclci. The level den- 
sity is less for the samarium nuclei indicating less single 
particle orbitals in the vicinity of the Fermi-energy. This 
is understood in terms of shell effects: the samarium nu- 
clei, with = 86 and 87, are closer to the A^ = 82 
neutron gap which means there are fewer orbitals below 
the neutron Fermi-energy. 

After establishing the level density as a function of 
excitation energy, we can explore various thermodynam- 
ical properties of the nucleus j|J^,^ . The experimental 
level density p{E) is proportional to the number of states 
accessible to the nuclear system, in this reaction, at ex- 
citation energy E. Nuclei are isolated systems with well 
defined energy, so the proper way to describe their sta- 
tistical properties is to use the micro-canonical ensemble. 
However, the canonical ensemble, permitting heat ex- 
change, and the grand-canonical ensemble, which in addi- 
tion allows particle exchange, are often used due to math- 
ematical difficulties with detailed calculations within the 
micro-canonical ensemble. In this work, thermodynamic 
properties of the ^''^'^''^Sm nuclei will be extracted and 
discussed using both the micro-canonical and the canon- 
ical ensemble. 



A. Micro-canonical ensemble 

When using the micro-canonical ensemble the partition 
function is simply the multiplicity of nuclear states uj{E), 
which experimentally corresponds to the level density of 
accessible states. Thus, the entropy is defined as: 



S{E) = kB \nuj{E), 



(8) 
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where u!{E) = p{E)/po and fcs is the Boltzmann con- 
stant. The normahzation constant po can be adjusted to 
fulfill the condition of the third law of thermodynamics; 
S —^ when T — > 0, T being the nuclear temperature. 
The ground-state band of even-even nuclei is assumed to 
have T = 0. Therefore, 5*0 is determined so that the en- 
tropy of the ground-state band in ^*^Sm is approximately 
zero. From Fig. | we see that p{E = 0) = 0.3 MeV-^ 
which determines po — 0.3 MeV~^. The entropy curve 
plotted on a linear scale is almost identical to the level 
density curve on a logarithmic scale as shown in Fig. |^. 
The difference in entropy between the Sm nuclei 

SS{E) = kg lnp(£;/^^Sm) - kg In p{E,^'^^Sm) (9) 

seems to be constant and equal ^ l.S/cs in the excitation 
energy region between E — 1.8 - 4.8 MeV. This has also 
been seen before in the well deformed nuclei (25) . As dis- 
cussed in detail in Ref. 6S{E) can be understood as 
the entropy of the unpaired neutron. 

The micro-canonical temperature in units of MeV is 
given by: 

The deduced temperature spectra exhibit pronounced 
bump structures, since the small bumps in the entropy 
curves are enhanced through the differentiation per- 
formed in Eq. (|l0|), as shown for i48a49gj^ 
per panel of Fig.l^ The bump structures are interpreted 
as the breaking of nucleon pairs. When particle pairs 
are broken new degrees of freedom open up i.e. a more 
than normal opening of new domains of the phase space, 
leading to an unusual increase in the entropy and thus 
the decrease in the temperature. We therefore interpret 
that the location of the break up of a nucleon pair hap- 
pens in the region where the micro-canonical temperature 
curve has a negative slope. From Fig. Qwe find that this 
takes place at ~ 1.9 MeV for ^^^Sm. This value can 
be compared to 2A„ and 2Ap, which are the expected 
cost of breaking up a neutron or proton pair. The pair- 
ing gap parameters Ap and A„ can be determined from 
empirical masses of a sequence of isotones or isotopes 
|6|. This gives 2Ap = 2.1 MeV and 2A„ = 2.2 MeV 
for ^''^Sm, which is close to 1.9 MeV found from Fig. |^. 
For -'^"'^Sm, the corresponding values are 2Ap = 2.7 MeV 
and 2A„ — 2.0 MeV while the first negative slope in 
the temperature curve of Fig. |^ is at about 1.75 MeV 
and the second negative slope is at about 2.5 MeV. The 
first negative slope in the temperature as a function of 
excitation energy in ^''^Dy, ^''^Dy, ^^^Yb, ^^^yb 0] and 
^^^Er 1^ also coincides roughly with 2A. While there is 
a deviation from 2A in the localization of the break up 
of the first pair oihaarticles in the case of -'^^^Er As 
discussed in Ref. ||] this can be due to the influence of 
structural effects in the nuclei, such as e.g. the Fermi-level 
position in the Nilsson single-particle scheme, variation 



in the density of single-particle orbitals, and two quasi- 
particle couphngs to collective degrees of freedom. 

The nuclear heat capacity deduced within the frame- 
work of the micro canonical ensemble can be determined 
from differentiating the temperature : 

= ' (11) 

and is shown in Fig. ^ for ^^^Sm and ^'^^Sm. The dou- 
ble differentiation of the entropy introduces strong fluc- 
tuations in the heat capacity and the decrease in the 
micro-canonical temperature lead to the spectacular fea- 
ture of negative heat capacity. The negative heat capac- 
ity has been seen before in ^^^Er and ^^^Er and discussed 
in Ref. ||] . It has also recently been observed experimen- 
tally in the critical region of nuclear fragmentation in Au 
quasi-projectile sources formed in Au-|-Au collisions [ p7| . 
There is some controversy on how to interpret negative 
heat capacities in small systems. Negative heat capacity 
is usually taken as a signal of a flrst order phase transi- 
tion suggesting that the process of breaking a nucleon 
pair around E ^ 2A appears like a phase transition of 
first order. Others claim that negative heat capacity in a 
small system like the nucleus is compatible with a second 
order phase transition [ p9| |. 



B. Canonical ensemble 



The transformation from the micro-canonical to the 
canonical ensemble is performed by the canonical parti- 
tion function 

oo 

z(r) = ^^ii;)e-^/^. (12) 

E=0 

The partition function is a Laplace transform of the mul- 
tiplicity of states i^{E) = p{E)/pQ. The experimental 
level density is only covering the excitation energy re- 
gion from zero to about Sn — ^ MeV. In the region of 
and above the neutron-separation energy Sn , the Fermi- 
gas model is believed to describe the nuclear properties. 
Therefore, the experimental level density is extrapolated 
to higher excitation energies with the Fermi-gas approxi- 
mation of Ref. [p3|[, (see the solid lines of Fig. |^) in order 
to evaluate Eq. (|12|) . The thermal average of the excita- 
tion energy in the canonical ensemble is 

oo 

{E{T)) = Z-^ ^ Eu{E)e-'^l'^. (13) 

By the Laplace transform in Eq. (jl^) much of the in- 
formation contained in the micro-canonical level density 
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becomes smeared out^. The lines in the upper panel of 
Fig. 1^ display the smooth variance of the canonical tem- 
perature as a function of the thermal average of the ex- 
citation energy. 

The nuclear heat capacity deduced within the frame- 
work of the canonical ensemble can be used as a "ther- 
mometer" for the quenching of pair correlations. It is 
the derivative of the thermal average of the excitation 
energy: 



Cv{T) = kB{l-T/Ty 



(15) 



Cv{T) 



d{E) 
' dT ' 



(14) 



The canonical heat capacities for ^^*Sm and ^"^^Sm are 
shown in Fig. ^ and exhibit an S-shape as a function of 
temperature indicating a second order phase transition 
0. The micro-canonical ensemble gives a detailed de- 
scription of the breaking of one, two, three,... nuclcon 
pairs as a function of excitation energy, while the canon- 
ical ensemble reveals the general average properties of 
this phase-like transition. The proper definition of ther- 
mally driven first- and second-order phase transitions in 
systems with few particles is a long standing problem, 
which will not be discussed in the present experimental 
work. 

The level density curve of ^*^Sm is more linear in a 
log-plot than the one for ^''^Sm (Fig. this leads to a 
more pronounced S-shape in the canonical heat capacity 
of the odd ^^^Sm nuclei compared to the even ^'^^Sm nu- 
clei, as shown in Fig. |^. At this point there is a striking 
difference between the samarium data presented here and 
earlier data on well deformed mid-shell rare earth nuclei. 
The opposite odd-even effect has been observed exper- 
imentally earlier between ^^^Dy and ^^^Dy, ^^^Yb and 

Yb d, and ^^^Er and ^^Tgr ||], where the even nu- 
clei show a more pronounced S-shape than the odd mass 
isotopes. Clearly there are different mechanisms at play 
here. In the region close to the TV = 82 close shell the 
odd nuclei are more deformed than the even nuclei. With 
deformation more single particle orbitals overlap giving 
more quasi-particle pairs the chance to interact, which 
again can be the reason for the sharper rise in the heat 
capacity curve giving a more pronounced S-shape. The 
differences in the S-shape of the heat capacity curves is 
an interesting observation which awaits further theoreti- 
cal efforts. 

The critical temperature for the quenching of pair cor- 
relations can be found by fitting the canonical heat ca- 
pacity with a constant temperature expression at low ex- 
citation energy 0. Assuming a constant temperature 
level density gives this expression for the canonical heat 
capacity: 



where r is the constant-temperature parameter identi- 
fied with the critical temperature. The dash-dotted lines 
of Fig. ^ describe Eq. ( p^ ) with asymptotes at r = Tc. 
The critical temperature for the quenching of pair cor- 
relations is found at = 0.45(5) MeV for ^''^Sm and 
Tc = 0.48(5) MeV for ^''^Sm, in the same range as for 
the deformed rare earth nuclei 0. The difference in Tc 
is within the error bars. 



V. DISCUSSION OF THE 7-STRENGTH 
FUNCTION 

There are several models developed for the 7-ray 
strength function fxL- We assume that the 7 decay in 
the continuum is dominated by dipole transitions. In 
Ref. j|] various theoretical models for El and Ml radia- 
tion were tested and compared to experimental data. In 
this work only the El and Ml models, which gave the 
best description of the data, will be used to model the 
^^^Sm and -'^^^Sm 7-strength functions. 

The model of Kadmenskii, Markushev, and Furman 
(KMF) is adopted to account for the El radiation: 



fEl{E.y) 



1 0.7aBl^|l(£;2-|-47^2^2^ 



Eei{E^ 



E? V 



(16) 



This model takes into account the energy and tem- 
perature dependence of the GEDR width and is often 
utilized to describe experimental data |jl^ . To be consis- 
tent with recent analysis performed for deformed nuclei 
the temperature parameter T is utilized here as a 
constant fit parameter. This assumption of a constant 
temperature is in accordance with the Axel-Brink hy- 
pothesis, Eq. (^, giving the same strength function for 
all excitation energies. The values for the giant electric 
dipole resonance parameters cr_Bi, and Eei for ^''^Sm 



are taken from |^l|], for ^*^Sm the values are unknown so 
the ones for ^''^Sm were used. 

The Ml radiation is described by a Lorentzian based 
on the existence of a giant magnetic dipole resonance 
(GMDR) 1,0: 



1 



STT^n^c^ {E^ - El,,Y + E^^T\,, ■ 



(17) 



where ctmi, Fmi and Emi are the giant magnetic 
dipole resonance (GMDR) parameters, which are taken 
from pi| . There are experiments which indicate such an 



'^The complete information is recovered when using a com- 
plex number for the temperature and thus defining the canon- 
ical partition function in the complex temperature plane. 
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Ml giant resonance due to spin-flip excitations in the 
nucleus |]3^ . 

The pygmy resonance is described with a similar 
Lorentzian function /py as Eq. (|l^), where the pygmy- 
resonance strength tTpy, width Fpy and centroid Epy 
have been fitted in order to adjust the total theoretical 
strength function 



/ = KifEi + L 



Ml J 



fp 



(18) 



to the experimental data. The resulting theoretical 7- 
strength functions are shown as solid lines in Fig. ^. The 
obtained values of the fitting parameters for the pygmy 
resonance are together with the normalization constant 
K given in Table The temperature parameter was first 
kept as a free parameter in the fit, giving T fa 0.5 MeV 
for ^■^^Sm, however the fit did not converge in the case of 
^^^Sm. Therefore, in order to be consistent we fixed the 
temperature parameter at 0.5 MeV in both cases and a 
more pronounced minimum was obtained. The cho- 
sen temperature of T = 0.5 MeV is in good agreement 
with the average temperature found in the micro canoni- 
cal ensemble as shown in Fig. ^ However, as it has been 
discussed in Q the temperature parameter in (16) may 
not correspond to the real nuclear temperature, due to 
possible uncertainties in the temperature dependence of 
the GEDR width. 

Figure ^ also shows the predicted individual contri- 
butions from the giant electric dipole resonance fsi, the 
giant magnetic dipole resonance /mi and the pygmy reso- 
nance fpy to the total 7-strength function. The strength 
function is generally dominated by El radiation. The 
Ml-strength function /mi is ~ 20% of the £'l-strength 
function f^i and the pygmy resonance contributes con- 
siderably around 3 MeV. 

The 7-strength functions of ^^^Sm and ^"'^Sm are com- 
pared with the 7-strength functions of the well deformed 
I6i,i62i3y i66467Er | and i7i472Yb g nuclei in 
Fig. 10. The 7-strength functions of ^''^Sm is less than for 



^'^^Sm. Some of the explanation can be that in the tran- 
sitional region from deformed to spherical shapes, one 
neutron can make a significant difference for the nuclear 
charge distribution. 

The pygmy resonance in the 7-strength functions was 
first explained in Ref. |Q by the enhancement of the 
El 7-strength function. This interpretation was adopted 
in Ref. M for the pygmy resonance found at ~ 3 MeV 
in I6i462j-)y g^^^ i7i,172yi-,^ According to systematic in- 
vestigations of the pygmy resonance parameters in the 
rare earth nuclei |Q it should appear around 1.7 MeV 
for ^''^'^^^Sm. However, the possibility that the pygmy 
resonance is of Ml character can not be excluded: at 
an excitation energy around 3 MeV, there is a concen- 
tration of orbital Ml strength in the weakly collective 
scissors mode [Q. The scissors mode was first observed 
in electron-scattering experiments , and is confirmed 
by the (7,7') reaction ||3^. Recent investigations of Dy 
and Yb isotopes M indicate that the observed bump 



can be interpreted as due to the scissors mode built on 
exited states. The present 7-strength functions show a 
bump around 2.5 MeV for ^^^Sm and around 2.0 MeV 
for "'^'^^Sm, as shown in Fig. The pygmy resonance in 
^^^Sm 7-strength function is less obvious. It might be 
that shell effects in ^'^^Sm hide the fine structure in the 
7-strength function. Also in the (n,7)-data for ^""^Sm in 
Fig. 1^ it is difficult to locate a 7-ray bump in the 3 
MeV region. 

In Ref. the total Ml strength, observed experimen- 
tally in (7,7') reactions, is found to be: ^i?(Ml) 1— 
0.43/i^ for ^''^Sm. Our data can be compared to this 
result using 



' ^ ^ ' £:py 327r2 



(19) 



For details see Ref. Hi. We obtain 3.5 for i^^^Sm and 
2.2 nj^ for ^''^Sm. The difference of a factor 10 between 
the results for ^**Sm can be understood by the fact that 
in Ref. [!§ only 7-rays between 2.7 MeV and 3.7 MeV 
are taken into account. They are therefore only summing 
up a part of the high energy tail of the resonance. In ad- 
dition weak 7-lines probably have escaped detection, as 
discussed in Ref. [ p8[ . 

The pygmy resonance parameters are shown in Fig. |ll| 
and we believe that the pygmy resonance in the samar- 
ium data is of the same nature as the resonance observed 
in the other rare earth nuclei. Looking at Fig. |l^ there 
seems to be a small odd-even effect in the strength of the 
pygmy resonance dpy. Clarification of the multipolarity 
of this pygmy resonance awaits new experimental data. 

For the first time a comparison of our data to the 
GEDR data ||l[ is performed, see Fig. [l|. The tail of 
the GEDR was until now not very well known and has 
been modeled in several ways. Our experimental data 
represent a unique guideline to theoretically describe the 
shape of this tail. 



VI. CONCLUSIONS 

The level density and the 7-strength function have 
been extracted experimentally for ^''^Sm and ^''^Sm. The 
level density of the samarium nuclei is smaller than for 
the mid-shell dysprosium, erbium and ytterbium nuclei 
investigated earlier. From the level densities thermo- 
dynamical quantities such as temperature and heat ca- 
pacity were found. Structures in the micro-canonical 
temperature are interpreted as the onset of new degrees 
of freedom by the breaking of Cooper pairs. The S-shape 
in the heat capacity curves, found within the canonical 
ensemble, indicates the pairing-phase transition, and a 
critical temperature for the quenching of pair correla- 
tions is found. The S-shape in the heat capacity curve 
is slightly more pronounced in the odd ^*^Sm nuclei, the 
opposite of what has been observed in the well deformed 
nuclei, where the even isotopes exhibit more pronounced 
S-shapes than the odd isotopes. 
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For the total 7-spectra of the ^^^Sm and ^^^Sm 
nuclei, from excitation energy at the neutron bind- 
ing energy, there is very good agreement between the 
neutron capture data from the ^^^Sm(n,7)^^*Sm and 
^^^Sm(n,7)^^^Sm reactions and the same spectrum cal- 
culated from the 7-strcngth functions and level densi- 
ties extracted from the present ^'**Sm('^He,a)^^^Sm and 
"8Sm(3He,3He')^'*^Sm data. 

The interpretation of the fine structure of the 7- 
strength functions is still an open question, since there 
are no experimental data to determine if the multipolar- 
ity of the pygmy resonance is of Ml or El nature. If it 
turns out to be Ml, the resonance at ~ 3 MeV in the 
7-strength function is likely to be interpreted as due to 
the scissors mode built on exited states. 

For the first time we compare our data to the (n,7) 
GEDR data and the ^'^^Sm 7-strength function seems to 
give a consistent extrapolation down to lower 7-energies. 
The tail of the GEDR was until now not very well known 
and has been modeled in several ways, and our experi- 
mental data will assist the theoretical effort to describe 
the shape of this tail. 
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TABLE I. Parameters used in the normalization of the level densities and 7-strength functions. The values of the different 
parameters are taken from ^ Ref. ||l[, ^ Ref. [§o) " Ref. @ and ^ Ref. H]. 



Nucleus 


P(B„) 
(10^MeV~i) 


B„ 
(MeV) 


a 

(MeV"i) 


D 

(eV) 


/ 


(r> 

(meV) 


"«Sm 
i*«Sm 


1.59(18) 
0.49(11) 


8.14 
5.87 


19.18^ 
19.85=^ 


5.7(5)^ 
100(20)1 


1.1 
1.6 


69(2)1 
45(2f 



TABLE IL Parameters obtained from the pygmy resonance in the 7-ray strength function. The fitting procedure is performed 



with the temperature as a fixed parameter of T = 0.5 MeV. 



Nucleus 






Tpy 


K 




(MeV) 


(mb) 


(MeV) 




i^«Sm 


1.99(8) 


0.08(3) 


2.5(2) 


0.61(6) 


"«Sm 


2.46(5) 


0.11(5) 


1.4(2) 


1.1(3) 
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FIG. 1. Gamma-ray spectra at 5 MeV of excitation energy for Sm (top) and Sm (bottom). The raw (left), unfolded 
(middle) and primary (right) 7-ray spectra are shown. 
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FIG. 2. The experimental level density for ^''*Sm (data points) is normalized between the arrows to known levels from 
spectroscopy [ p^ at low excitation energy (histogram upper panel) and to the level density calculated from neutron-resonance 
spacings [JL9| at the neutron separations energy (square data point in lower panel). The solid line is a Fermi-gas approximation. 
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FIG. 3. The level densities for ^**Sm and ^^^Sm (data points). The solid lines are Fermi gas approximations for the level 
densities. The filled squares are the level densities at the neutron binding energy calculated from neutron-resonance spacings 
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FIG. 4. The 7-strength functions for "*Sm (circles) and "^Sm (triangles). The error bars are due to statistical uncertainties, 
only. 
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FIG. 5. The total 7-spectrum for ^'**Sm and ^''^Sm measured from excitation energy equal to the neutron binding energy. The 
data points with error bars are taken from ^'*'^Sm(n,7)^'**Sm and ^'**Sm(n,7)^'*®Sm experiments The sohd hne is calculated 
from the 7-strength functions and level densities extracted from the present ^'**Sm(^He,Q)^^^Sm and ^''*Sm("^He,^He')^**Sm 
reaction. The calculation is performed by averaging over 100 keV intervals. 
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FIG. 6. The level density functions for "**'i''^Sm (filled triangles), "^-i^^Dy [| (open triangles), i^^'ie^Er |§ (filled circles) 
and i7i,i72Yb g (stars). 
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FIG. 7. Temperature and heat capacity for ^''^Sm and ^*^Sm derived within the micro canonical ensemble (filled circles) and 
canonical ensemble (solid lines). 
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FIG. 8. The heat capacity for "**Sm and Sm as a function of temperature derived within the canonical ensemble. The 
dashed dotted Hne describes the constant temperature estimate, Eq. (^5|), where t is recognized as the critical temperature and 
marked with the vertical line. 
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FIG. 9. The experimental 7-ray strength functions (data points) of ^"^^Sm (left) and ^'*^Sm (right). The solid line is the fit 
to the data by the theoretical model. The dashed lines are the respective contributions of the GEDR fsi, the GMDR /mi, 
and the pygmy resonance fpy to the total theoretical strength function (solid line) . 
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FIG. 11. Systematics of the pygmy resonance parameters for odd (filled squares) and even (open squares) nuclei as a function 
of neutron number N. The resonance energy Epy, the width Fpy and the cross sections apy with error bars are shown in the 
upper, middle and lower panels, respectively. 
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FIG. 12. The experimental 7-strength function for '^^^Sm {fsi + /mi) from the present (''He,a)-reaction (circles) plotted 
together with data obtained from photo-absorption cross section Ref. M (triangles). 
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